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Abstract 

An explicit formula for the chaotic representation of the powers of increments, (Xt+t — Xt ) n , 
of a Levy process is presented. There are two different chaos expansions of a square integrable 
functional of a Levy process: one with respect to the compensated Poisson random measure and 
the other with respect to the orthogonal compensated powers of the jumps of the Levy process. 
Computationally explicit formulae for both of these chaos expansions of (Xt+t ~ Xt ) n are given 
in this paper. Simulation results verify that the representation is satisfactory. The CRP of a 
number of financial derivatives can be found by expressing them in terms of (Xt+t — Xt ) n using 
Taylor's expansion. 

MSC: 60J30; 60H05 

Keywords: Chaotic representation property; Levy process; Power jump process; Poisson random 
measure; Martingale representation. 



STATISTICS SECTION TECHNICAL REPORT TR-07-04 



2 



1 Introduction 



The chaotic representation of a square integrable functional of a Le vy process i s an 



Sole et al. 



expan- 



(|2006l N l for a 



sion via its expectation plus a sum of iterated stochastic integrals, see 
recent review of such representations. There are two different types of chaos expansions: 



Ito 



(jl956l ) proved a Chaotic Representation Property (CRP) for any square integrable functional for 

e integrals with respect to a 



Nualart & Schoutensl (|2000h 



a general Levy process. This representation is written using multip 
two-parameter random measure associated with the Levy process, 
proved the existence of a new version of the CRP, which states that every square integrable Levy 
functional can be represented as its expectation plus an infinite sum of stochastic integrals with 
respect to th e ortho gonal ized compensated power jump processes of the underlying Levy process. 



Benth et al. 



(|2003h and 



Sole et al. 



(2006) derived the relationships between these two repre- 
sentations. However, these representations are computationally intractable. For the powers of 
increments, (^+t„ — X to ) n , of a Levy process, we instead derive computationally explicit formulae 
for the integrands of these two chaotic expansions. Hence we have all the results necessary to 
construct arbitrarily accurate computational formulae for the Levy functionals themselves. 



Pow er jump processes are important in mathematical finance. iBarndorff-Nielsen &: Shephard 
(|2006l ) performed hypothesis tests on exchange data under the null of no jumps and found that the 
tests were rejected frequently. In fact, at intraday scales, prices move essentially by jumps and 
even at the scale of months, the discontinuous behavior cannot be ignored in general. Only after 
coarse-graining their behavior over longer time scales do we obtain something similar to Brownian 
motion. Jumps can be understood both in te rms of a Poisson random m easure, or equivalently, 
by using the Power jump processes. Note that lNualart fe Schoutensl (|2000l . Proposition 2) proved 



that all square integrable random variables, adapted to the filtration generated by the Levy process 
denoted by X, can be represented as a linear combination of powers of increments of X, see Section 
12. II below. In fact, for any square integrable random variable, F, with derivatives of all order, we 
can apply Taylor's Theorem to express F in terms of a polynomial of powers of increments of X. 
Thus, the CRP of a number of financial derivatives can be found using this method, as is discussed 
further in Section [5j 

The deriv ation of an explicit formula f o r the CRP has been the focus of co nside rable study, see 



Leon et al. 



(120021 ) . lLokkal ([200J) and 



Eddahbi et al. 



for ex ample iNualart Sz Schoutensl (|200ll ) , 
(|2005l ). All the explicit formulae for general Levy functionals derived in these papers use the 



Malliavin type derivatives to derive explicit representations of stochastic processes for applications 
in finance. The derivative operator D is, in all of these cases, defined by its action on the chaos 
expansions. In other words, the explicit chaos expansion must in fact be known before D can 
be applied to find the explicit form of the pred i ctable or chaotic representation, thus yielding a 
circular specification. For example. iLeon et al. I (J2002J, Definition 1.7) defined the derivative of F 
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in the /-direction by 
,(0 



n=l ii,...,i n k=l 



(*)(•) 



and 



Leikkal (|2004l . Section 3) defined the derivative operator by 

oo 

D t]Z F = Y,nI n -i (fn{; t,z)) 



n=l 



where 



In (fn) 



fn —,t n , Z\ j .. . , Z n ) d(/i - 7r) 



'[0,T] n xL 

(Please refer to the corresponding papers for notation). Note that both of these definitions require 
the knowledge of the functions {fi lt ...,i n Ys or /„ (ii, t n , Zi, z n ) : s, which are the integrands of 
the chaos expansion of F. 



Jamshidian 



(|2005h extended the CRP in 



Nualart &: Schouten; 



(2000) to a large class of semi- 



martingales and derived the explicit representation of the power of a Levy process with respect to 
the corresponding non-compensated power jump processes, which is discussed futher after Theo- 
rem H] in this paper. Note that Le vy processes are included in the class of semimartingales, see 

, Corollary 2.3.21, p. 92). Our formula gives the explicit repre- 



Kannan Lakshmikantham (2001 



sentation with respect to the orthogonalized compensated power jump processes. Our result is 



therefore complementary to Jamshidian's formula, si nce our explicit formu 
respect to the orthogonalized processes, as defined by 



Nualart Schoutens 



a give s the CRP with 
(|200oI l 



In practical applications, it is often convenient to truncate the representation given by the 
PRP. The truncated representation of a stochastic process would yield a practically implementable 
approximation to the stochastic process. This approximation would be used for simulating the 
process, or with a finite number of traded higher order options, providing pricing formulae. The 
truncation would be chosen with minimal variance constraint. The advantage of expressing the 
sum in terms of stochastic integrals with respect to the orthogonalized processes is that the error 
terms omitted will be uncorrelated with the terms remained in the approximation. Jamshidian's 
result holds for general semimartingales (a larger class than ours) but our formula is designed for 
those with compensators equal to a constant times t only (which is satisfied by all Levy processes). 
Our results can be easily extended to semimartingales when the form of the compensators is 
known. 



The rest of the paper is arranged as follow: Section [5] gives the background information about 
the CRP for Levy processes. We give the explicit formulae for the CRP for (X t +t — X to ) n of 
a Levy process in terms of power jump processes in Section [3] and in terms of Poisson random 
measure in Section [U We show that in the Levy case, our formula complements Jamshidian's 
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formula. Section [5] gives the representation of a common kind of Levy functionals with the use of 
Taylor's Theorem. Simulation results for the explicit formulae are given in Section [6l In Section 
some concluding remarks are provided. Proofs and plots are included as appendices at the end. 

2 Background 

2.1 Levy processes and their properties 



We give a brief account of Levy processes and refer the reader to the work by ISatol ()1999l ) for 
a more detailed account. A real- valued cadlag stochastic process X = {X t ,t > 0} defined in a 
complete probability space (fi, T, P) on R d is called Levy process if X has stationary and indepen- 
dent increments with Xq = 0, where T is the filtration generated by X : Tt = o~ {Xt , < s < t} . 
Denote the left limit process by Xt- = lim s ^t,s<t X s , t > 0, and the jump size at time t by 
AX t = X t -X t _. 

A Levy process is fully specified by its characteristic function. Let <f>x x (u) be the characteristic 
function of the Levy process at t = 1, X\, that is, 4>x 1 ( u ) = E [e 1 ^ 1 ] . The characteristic function 
of X t is the n given by ((f>Xi since the distribution of a Levy process is infinitely divisible, see 



Satd (|1999l . chapter 2) . The cumulant characteristic function ip (u) = log (px 1 (u) is often called 



the characteristic exponent, which satisfies the Levy-Khintchine formula: 



1 



+oo 



ip(u) = i-yu - -a u + / (exp (iux) - 1 - iuxl { \ x \ <1} ) v (dx) , (1) 

^ J -oo 

where 7 G R, a 2 > and 1^ is a measure on R\ {0} with v ({0}) = and 

/+00 
(l A x 2 ) v (dx) < 00. 
-00 



In general, a Levy process consists of three independent components: a linear deterministic com- 
ponent, a Brownian component and a pure jump component. The Levy measure v (dx) dictates 
the jump process: jumps of sizes in set A occur according to a Poisson process with intensity 
parameter J A v (dx). To model a generic Levy process, only 7, a and a form for v (dx) need to 
be specified. 

In the rest of the paper, we assume that all Levy measures concerned satisfy, for some e > 
and A > 0, 

exp (A |x|) v (dx) < 00. (2) 



(-e,eT 

This condition implies that for i > 2, J_ 00 |x|*z;(dx) < 00, and that the characteristic function 
-E [exp (\uXt)} is analytic in a neighborhood of 0. 

Denote the i-th power jump process by = ^ 0<S<( (AA S ) 1 , i l> 2, and for completeness let 
= Xt. In general, it is not true that Xt = X^o<s<t ^X s ; this holds only in the bounded vari- 
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ation case, with a 2 = 0. By definition, the quadratic variation of X t , [X,X] t = ^ 0<s<t (AA" s 



X t when a 2 = 0. The power jump processes are also Levy processes a nd jump at the same 



Nualart fc Schoutens 



time a s Xt, but with jump sizes equal to the i-th powers of those of Xt, see 

(hood). 

Clearly -EpQ] = SpQ 1 ^] = m\t, where m\ < oo is a constant and by Protterl ([20041 p. 32), we 
have 



E[xf ] \ = E[ { Ax sY} = t I x i u(dx) = nut < oo, for i > 2, 



0<s<t 



(3) 



thus defining mi. iNualart Schoutens! (|2000l ) introduced the compensated power jump process 
(or Teugels martingale) of order i, ^Y®\ , defined by 



Y t (i) = X y t l> - E[X y t l >] = Xf 1 - mit for i = 1,2, 3, .... 



-CO 



(0 



(4) 



Nualart Schouten; 



li) 

Y t is constructed to have a zero mean. It was shown by 
that there exist constants an,ai 2, ai j_i such that the processes defined by 



(|2000l . Section 2) 



(5) 



for i > 1 are a set of pairwise strongly orthogonal martin gales, and this implies that for i ^ j, 
the process i s a m artingale, see 



Leon et al. 



([2002]). For convenience, we define a^i = 1. 



Nualart &: Schoutens! (|2000l ) proved that this strong orthogonality is equivalent to the existence of 
an orthogonal family of polynomials with respect to the measure 



d?7 (x) = O" 2 d5o (x) + x 2 u(dx) 



where 5q (x) = 1 when x = and zero otherwise, that is, the polynomials p n defined by 



Pn (x) = ^2 a n jx j 1 



are orthogonal with respect to the measure 77: 



Pn (x) p m (x) dr) (x) =0, n^m. 



We now state some key re lated results in the representation of stochastic processes given in 
Nualart fc Schoutens! (|200oh . 



Denseness of polynomials ([Nualart &; Schoutens! ([20001 . Proposition 2)): Let V = {X^ (X t2 — 
X tl ) k * ■ ■ ■ (X tn - X tn _J fc « : n > 0, < h < t 2 < ■ ■ ■ < t n , h, ...,k n > 1} be a family of 
stochastic processes. Then we have that V is a total family in L 2 (Q.,J r T,P) , that is, the 
linear subspace spanned by V is dense in L 2 (Q,,J-t,P); each element in L 2 (f2,.F) can be 
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represented as a linear combination of elements in V. 

• Chaotic Representation Property (CRP): Every random variable F in L 2 (f2,.F) has a repre- 
sentation of the form 

F = E(F) + Y^ E / / -/ /(fe,... A )(ti,.--,*i)dfff ) ---djHr^djff4* l) > (6) 
, i ; : ,. -i 7 " 7o 7 ° 

where the /(i lv .. i i i )'s are functions in L 2 (R^). This result means that every random vari- 
able in L 2 (p,,J r T,P) can be expressed as its expectation plus an infinite sum of zero mean 
stochastic integrals with respect to the orthogonalized compensated power jump processes 
of the underlying Levy process. Note that this representation does not explicitly allow for 
calculation of the integrands. 

• Predictable Representation Property (PRP) : From the CRP stated above, we note that every 
random variable F in L 2 (Q,^, P) has a representation of the form 

°° POO 

F = E[F] + Y J #dtf«, (7) 
t =i 7o 

(i) 

where <j) s s are predictable, that is, they are ^--measurable. 



2.2 Jamshidian's notation 



In 



Jamshidianl (|2005l ). which extends the CRP to semima rtingales, the power jump pr ocesses 



and compensators were denoted and defined diffe r ently from iNualart fc Schoutend (|2000l ). The 
power jump processes were defined in I Jamshidianl (|2005l ) by 



[X] 



(2) 



[X c ] t + (AX S ) 2 and [Xp = (& X s) n for n = 3, 4, 5, 



(8) 



s<t 



s<t 



(2) 

where [X c ] t = [X]^ is the continuous finite- variation (not martingale) part of [X]\ . Note that 
Jamshidian suppressed the time index t, but we add it here for clarification. Jamshidian denoted 
the compensator of by {X)^ . The compensator, (X)^ , is the predictable right-continuous 



finite variation process such that [X][ n ^ — (X)f l> is a uniformly integrable martingale. The 



compensated power jump process, denoted by X, 



(n) 



(«) 



[X] 



(«) 



(X)^ forn 



is thus defined by 
2,3,4,.... 



(9) 



For Levy processes, the compensators have the form mrf, where m\t = E[Xt] and mrf = 
t f^° Qo x l i'(dx) for % = 2,3,4,.... However, for semimartingales, the general form of the compen- 
sators is not known. 
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3 The chaotic representation with respect to power jump processes 

In this section we firstly derive the explicit formulae for the CRP when F in is the power 
increment of a pure jump Levy process and extend it to a general Levy process. 

3.1 Pure jump case 

Let us first outline the form of the representation to introduce the reader to the flavour of results 
in this paper. Suppose to > an d let Gt be a pure jump Levy process with no Brownian part 
(that is, a 2 = 0), be its i-th power jump process and G^ be its i-th compensated power jump 
process. Based on the structure of the expressions for (Gt+t — Gt ) 3 and (Gt+t — Gt ) 4 calculated 
using (|A.4|) - (|A.6|) . we desire to derive a general formula for (Gt+t — Gt ) k , fe = 1, 2, 3,..., as this 
forms a starting point for the representation of Xf. We notice that the number of stochastic 
integrals in eac h of the above representatio n is less than the possible full representation specified 
in the CRP by 



Nualart fc Schoutensl (|2000h : 



rt-\-t{) ft\_— ptj — i — 

(x t+t0 -x t0 ) k = f^(t,t ) + Y, E / / ••/ a'-'-'i 

1 f - . \ J to to J to 



3=1 

e{i,...,k} j 
xd yfe)... d yfe) d y^) 5 

where the i ,^'s are deterministic functions in 1? [R+j an( i are defined in Q. In 

(G t+t0 - G t0 )\ we have J^ J^dG^dG^, // +to dG« and J^dG^, that we shall repre- 
sent via the list {(1, 1) , (1) , (2)} . We can do an equivalent representation of (Gt+t — Gt ) 3 and 
{Gt+t — Gt ) A to get the following two lists: 

{(1,1,1), (1,1), (1,2), (2,1), (1), (2), (3)}. 

{(1,1, 1,1), (1,1,1), (1,1,2), (1,2,1), (2,1,1), 

(1, 1) , (1, 2) , (2, 1) , (2, 2) , (1, 3) , (3, 1) , (1) , (2) , (3) , (4)} . 

In general, the list of the orders of the compensated power jump processes of the stochastic integrals 
in (Gt+t — Gt ) k depends on the collection of numbers 



Ik = < 



(ii,i 2 ,...,ij) | j E {1,2, ...,&}, i p G {1,2, ...,k} and ^ i p < k } . (10) 

P =i 



This construction is explained in the beginning of the proof of Theorem [T] (Appendix [C]) using 
induction. A typical element ij) in therefore indexes a multiple stochastic integral 

j-times repeated with respect to the power jump processes with powers and indexed 
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t\, t%, tj. That is, (ii,i2) indexes the integral 

"t + t Q ft\— ftj-l ,. s ,■ \ f. s 

/ •••/ dG^-dO^dG^. 

to ■/ to J to 

Next we consider the terms in the representation not involving any stochastic integrals. That is, 
in (Gt+t — Gt ) 2 , m\t 2 + rri2t is considered; in (Gt+t ~ Gt ) 3 , mft 3 + 3mim2t 2 + m^t is considered, 
and in (Gt+t — Gt ) 4 , m\t 4 + 6m 2 m,2t 3 + (4mim3 + 3ml) t 2 + m&t is considered. We use (|A.4|) - 
()A.6p . given in Appendix |Aj to derive the representation. This time the representation can be 
simplified a great deal since we are not considering any stochastic integrals. Denote the terms 
which do not contain any stochastic integrals in (Gt+t ~ Gt ) k by c[ k \ 

Proposition 1 cj r) = for all r, cf ] = 1, G t (1) = nut, and for k = 2,3,4, 

C ^ = E ( •) ™i tC t j) - E (*) mi f h dCt 3) + m k t. (11) 
j—l V J / J=l 



Proof. The results for cj r) and C t (0) are trivial. For fc = 1, (G t+to - G to ) = J^dG^ + m\t 
and hence G^ = m\t. For k > 2, the terms in ()A.4p are equal to zero since Gt has no Brownian 
part (a 2 = 0). The first term in (|A.5|) contains a stochastic integral and hence from the second 
term of (|A.5|) and (|A.6[) . we have 

c * k) = E ( •) ^ (* + *>) ~ i} - E ( ■) m > r° dC tS + m k t. 

j=l \j ' j=i ^ ' J *o 

Putting u = t\ — to in the second term, we have 

c ( fc) = ^ A\ m . (t + k) c m _ ^ A\ m . /" (u + ^ dC ( fc _,) + mfct 

= E ( •) - E ( •) ^ jf f i + 

(fc) 

Note that G t is independent of to- D 



Thus, given Proposition [TJ G t can be expressed in terms of m/s for any given k and easily 
coded. We will show in the followings that in the calculation of (Gt+t — Gt ) k , all the G t s, j = 
0,1,..., A; are required. In fact the coefficients of the stochastic integrals in the representation 
depend only on C\ s, j = 0, 1, k, as stated in Theorem [JJ below. 
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(k) 

The next proposition gives the representation for C t in a non-recursive form. Let 



C-k = < 



(il,i2, k) \l G {1, 2, fc} , G {1, 2, A;} , it > i 2 > ■ ■ ■ > i\ and ^ i q = k } . (12) 

9=1 



The number of distinct values in a tuple (fik = (i^,^, in Ck is less than or equal to I. 

When it is less than I, it means some of the value(s) in the tuple are repeated. Let the number of 
times r G {1, 2, 3, .., k} appears in the tuple 4>k = if^ J be pt k ■ 

Proposition 2 

t' (13) 



, I .(k) .Ik) .lk)\^ r 



n 



where , are the elements of <pk, pt k, s are defined above and (i[ ,4 , )! is the 



multinomial coefficient: fij^,*^, zj^ J ! = ( fc)| j , 



. V fe) i-i (fe) ' 



Proof. Proof is included in Appendix [Bj □ 

Let nW A| .„ A)>t be the coefficient of jj" • • • J^^Gf > • • • d^dG^ in (G t+fo - G to )< 
We then have the following result. 



Proposition 3 



3 

n (ii,i2,...,ij).t = fa' 1 ' 2 ' ") !C t where n = fc "EV ( 14 ) 



I,.. .,<,]/ 

p=l 



Proof. The proof of Proposition [3] is contained in the proof of Theorem [TJ □ 

For example, say we want to find the coefficient of J t * + *° dG^dG^ in {Gt+t — Gj ) 4 , that 
is, we want to find ID* ^ v To derive this coefficient, we first note that n = 2 and so n[i\) t = 
i!i!2! ^ = 12 (^t + mf t 2 ) , which can be easily verified by calculating (Gt+t — Gt ) 4 using 
(|A.4j) - (|A.6|l . Now we put the above results together to get a general formula for {Gt+t Q — Gt ) k . 

Theorem 1 Let Gt be a Levy process with no Brownian part satisfying condition ([2]). Then 
the power of its increment can be expressed by: 

(G t+t0 - G to f = KW,t,t + C[ k \ (15) 
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where Ik is defined in (|10p , Yvq t is denned in Proposition [3j the C\ are constants defined in 
Proposition [2] and <S(i ll i 2 ,...,i ),t,io i s defined as the integral: 

J to J to J to 

Proof. Proof is included in Appendix [Cj □ 



To derive the explicit formula for the power of increment of a Levy process, {X t+to — X to ) n , with 
respect to orthogonalized compensated power jump processes, we need the following proposition. 

Proposition 4 The n-th compensated power jump processes, Y^ n \ of a general Levy processes 
satisfying condition ([2]), can be expressed in terms of the orthogonalized compensated power 
jump processes, for % = 1, 2, ...,n, by 

n-l 

fc=i 

where b n ^ denotes the sum of the set M. n,k , which is defined by 

M n ' k = a iu i 2 a i2 ,i 3 ■ ■■a ij _ lt i ] : h = n,ij = k, i p > i q if p < q, i p G N for all p| , 

and M n ' n = {1} . 

Proof. Proof is included in Appendix [Dj □ 



Theorem 2 Let Gt be a Levy process with no Brownian part satisfying condition ([2]). Then 
the power of its increment in terms of stochastic integrals with respect to the orthogonal 
martingales, H, is given by the following equation: 

(G t+t0 - G t0 ) k = £ U^X+C^, (16) 

where is defined in (jlOp , LT^ t is defined in Proposition [3j c[ k ^ is defined in Proposition [2] and 
sj-. ■ s . . is defined as the integral: 

S (H) 
(ii,t2,...,ij),t,to 

'i *3 /•*+<<) /-ti- rtj-i- , . , , 

= E--- E E^-WA*/ / ■■■/ d^-.-d^d^, 

6 nj fc is defined in Proposition [H 
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Proof. From Proposition 01 we have 



rt+to rti- 



to Jt 

t+to rti- 



tj-i- 



dG, (il) • • • dG): j " i; dG 



to 



t2 



•3-1- 



10 J to J to 

11 *3 



^ b ilikl H ( t 



(ki) 

3 



fel=l 



•••d 



(*J-i) 

t2 



fcj_l=l 



kjV 7 ' 



ti 



k\=l kj — i=l kj=l 

Hence, by using Theorem [H we finish the proof. 



t+t fti- 



■3-1- 



r(*i) 



dHr j ■■■dH^ j - l) dH^ kj) . 



□ 



Corollary 1 By TheoremUl 



(Gt+to - Gt ) m (G t +t - G to ] 



y u^\sf\ t +c\ m) I I V nf\sF\ t +c t 

/ J Vm,t Vm,t,tQ I I I / j V n ,t V n ,t,to 1 




Jn) AH) 



y nS m+n l5f } M +c t 



n £-£71 

(m+n) 



Hence, we can find out how to express the product of two iterative stochastic integrals of orders 
m and n as a weighted sum of iterative stochastic integrals of order m + n, m + n — 1,..., 2, 1. 

Note in Theorems Q] and El the integrands of the stochastic integrals do not involve to nor 

(p) 

any of the integrating variables t%, t2, tj. They are completely characterized by C t s, where 
p = 0,1, k. Hence to find the chaotic representation of (Gt+t — Gt ) k , we only need to know 
the moments of Gt, m\t = E [X t ] and m p = x v v (dx) for p = 2, k. This result is intuitive 
as (Gt+t — Gt ) is a stationary process. 

3.2 General case 

Next we want to derive the formula for the power of the increments of Levy processes with 
a ^ 0. Recall X = {X t ,t > 0} denotes a general Levy process, X t denotes its i-th power 
jump process and Y t denotes its i-th compensated power jump process as defined in fl4J). We 



define A\ (X t+to , X to ; k) and A 2 (X t+to , X to ; k) such that (X, 



t+k, 



X 



to i 



At (X t+t0 ,X t0 ;k) + 



A2 (Xt+t ' -^to i ^) ' wbere A\ (Xt+to : -^to i ^) comprises all the terms not containing a in (X 



t+t 



X, 



to J 



STATISTICS SECTION TECHNICAL REPORT TR-07-04 



12 



From (|Xij) - ([AT6|) . it may be noted: 

'to 



(X t+t0 - X to f = — k (k - 1) [ (X t+t0 - X t0 Y" t 



pt+to \ 
jf (s-t )d(X s -X t0 ) k - 2 ) 

k /k\ r t+t ° 

+ E / A 2 (X s -,X t0 ;k-j)dY s (» 

j = l \J / J to 
k ~ 1 fk\ 

+ X, ( • ) m i (* + *°) ^ 2 (-Xt+to. - j) 

j=i w 

^ /V\ ft+U) 

- E • H / sd[i 2 (I s ,I t0 ;fc- i)] + (X t+t0 , X t0 ; fc) . (17) 
Proposition 5 For any Levy process Xt satisfying condition ([2]), 

(X t+t0 - X t0 ) k = A 1 (X t+t0 ,X t0 ;k) + (k _2n)\ ri¥ a2nAl " 2n) 

n=l 

Proof. The proof uses the same techniques as in the proof of Theorem[TJ Note that A\ (X t +t j Xt ; p) 
where p = 1, 2, k, are given by Theorem [1] □ 

Proposition Ogives the formula of (Xt+t ~ Xt ) k in terms of a summation of A±, where [k/2\ + 1 
calculations of Ai are needed. The next theorem gives the formula in an alternative form which 
requires A± to be computed once only. 

(k) 9 (k) 

Definition 1 Let C\ J. be the terms obtained by replacing rri2 with vri2 + a in C\ (Proposition 
[2]) and IT. . . w be the terms obtained by replacing C?^ with ci*L in IT^ . . w 
(Proposition [3]). We then note the following theorem. 



Theorem 3 For any Levy process X% with cr 2 7^ and satisfying condition ([2]), the 
representation of (X t +t ~ Xt ) n is given by Theorem [1] with m 2 replaced by (m,2 + o" 2 ) , i.e. 



(X t+to - X to ) n - ^ U el,t,a S e n ,t,t + C t?cr > 



where X n is defined in (fTO]) . Ilg" t(T and C* t ™ are defined above and S'^ ia j 4 tQ is defined to be 
the integral: 



to ^ to -' to 
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Proof. We define a new class of power jump processes by 



X^=X^ + a\ 

X® = X[ j) for j = 1 and j = 3, 4, 5, . 



(18) 



We also define a new class of compensators 



m 2 t = (m 2 + a 2 ) t, 

rhjt = rrijt for j = 1 and j = 3, 4, 5, .... 
Hence, by definition, the compensated power jump processes, 

y W = xf ] - rhit = Xp - rmt 



= Y t for all i > 1. Therefore, the representation of (Xt+t — Xt ) in terms of the stochastic 

(i) . (i) ~(i) 

integrals with respect to Y t is the same no matt er we start from using X^: o r X^ . To calculate 
the expression using xf\ we use equation (2) in Nualart &: Schoutens ( 2000 ): 



(X t+to - X to ) k 



Ell) I ° (Xs- - X t0 ) k -i dX^ 

j — l \J / Jto 

+ ^k(k- 1) ((X t+t0 - x t0 ) k - 2 1 - J 8 d (X s+t0 - x U) ) k - 2 ^ 
£ (j) J to - Xt ) k ~ j axp + a -k{k- 1) j( (x, 

JL /h\ ft+to , -.2 /-t+to 

EQ/ (*.--Xh)*" J d*. 0) + Y*(*-i)j[ (x„_-x« (1 ) t - 2 d, 1 

g G) £ + '° <x - - dxi " + 6) f <*- - **' ) " d ( * 2 



(■-•+to)- ~ X to) k 2 ds 



0- 



By (|lHp . we have 



(x Wo - x t() ) fc = £ (*) (X s _ - x to f-i dip). 

Using exactly the same calculation as the one leading to (|A~4"] ) -(|AT6 ]) . we have 

k /k\ f t+t ° fe_1 / A 

(X t+t0 - X t() ) k = J2 (j) j tQ (Xs- - X t0 ) k ~i 6yW + J2 m," (* + *>) (X t+ to - X 
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This is exactly the equation (|A.5|) - (|A.6j) we based on in the derivation of Theorem [H except that 
nij is replaced by rhj. Hence we now have a simple formula for the representation of (Xt+t — Xt ) k 

(i) 

in terms of the stochastic integrals with respect to Y t by replacing rrij with rhj in the formula 
given by Theorem [TJ In other words, we have 

{X t+to - X t() ) n = ^2 u el,t,a S e n ,t,t + c £> 

(n) (n) 

where ' t and C^J are defined above. Note that this representation does not depend on the 

(i) 

power jump processes directly since it is in terms of the compensated power jump processes, Y t . 
So it does not matter if we change the definition of the power jump processes, as long as we change 
the compensators accordingly, we will get the same compensated power jump processes. □ 



Theorem 4 For any Levy process Xt with a 2 ^ and satisfying condition ([2]), the 
representation of {X t+ta — X to ) n is given by Theorem [2] with 771,2 replaced with (771-2 + & 2 ) > i-e. 



(v y \n _ ST^ ri{ n ) q'( H ) , r {n) 



where l n is defined in (I10D . T& 1 \ and C} n j are defined above and S' r ■ . w . 

" 4 »' a„,t,a t,cr (ii,i2,...,lj),t,to 

the integral: 



is defined to be 



S >(H) 

(ii,i2,...,ij),t,to 



ki=l fcj_i=l fcj=l 

6 nj fc is defined in Proposition [U 

Proof. It follows directly from Theorems [2] and [3l 



to J to 



to 



□ 



Remark As noted in Section 



Jamshidianl (|2005l ) derived an explicit formula for the chaotic 



representation of X* in terms of the non-compensated power jump processes, X\ J \ when X t 
is a semimartingale. Our explicit formula gives the representation in terms of orthogonalized 
compensated power jump processes, H^. In the following, we show that in the Levy case, our 
formula complements Jamshidian's one. We note the notation used by Jamshidian in Section 
12.21 If X is a Levy process, we can see that [X c } t = [X\ c t = a 2 t (where the superscript c stands 
for continuous part of the process) and hence [AT]^ = a 2 t + ^2 s<t (AA S ) 2 . With Jamshidian's 
notation, the a 2 is implicitly included in the [A^]j 2 \ 



Jamshidian 



(|2005l ) defined C = C* D C*, where C* is the set of semimartingales of finite moments 



with continuous compensators adapted to a Brownian filtration, and is the set of processes with 
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exponentially decreas ing law. Jamshidi an generalized the CRP from Levy processes to the set C. 
In proposition 8.2 of I Jamshidianl (|2005l ). an explicit formula for the chaotic representation with 
respect to the non-compensated power jump processes for the semimartingales in C when to = 
was derived. iJamshidianl (|2005l ) defined the power jump processes using the power brackets, see 



([8]) and Q. The multi-indices were denoted by I = ...,i p ) € N p , where N is the set of natural 
numbers, and for integers 1 < p < n, 



= i 1 = (»l,-i*p) G N p : ii + ■•• + i p = n} , p,neN. 



(19) 



k n 



Note that from (|1U|) . 2^ = |J |J N^. Proposition 8.2 of IJamshidianl (|2005l ) states that, for a 

n=l p=l 

semimartingale Xt with Xq = 0, we have, for all n £ N 



x " = E E 



k\---i P \ Jo Jo 



t ft!- 



tp—l — 



d[X]^...d[X^-''d[X]^>. 



(%>-!) a r vi(*j>) 



(20) 



Since IJamshidianl (|2005l ) only considered non-compensated processes, we substitute all the in 
IT]) by zeros (since the compensators in the Levy case are mji), which makes C$ = for all 



k t^O. Son 



(*) 



is non-zero only when X]p=i V = ^> as defined in (|19p . Hence in the Levy 



case, Theorem [3] reduces to (f20|) . 



Corollary 2 T/ie expectation of (Xt+t — Xt ) k is given by Ci" , which can be obtained by replac- 
ing ?7i2 with m2 + °" 2 , given by equation I113\) . 

Proof. As the expectations of all the stochastic integrals are zero, this follows directly from 
Theorem El □ 



Corollary 3 The expectation of 

7772 + o~ 2 and mi with in c[ k \ given by Proposition^ 



(i) 



can be obtained by replacing 777,2 with 



From Corollary [3 E [X?] can be obtained by replacing 7772 with 7772 + o 2 in C\ . Since 
= Xt — m\t and 



•/ 



by putting mi = in (|21j) . we can conclude that the expectation of 
by replacing 7772 with 7772 + <r 2 and 777 1 with in c[ k \ 



I'M? 



(21) 



can be obtained 



In the next section, we extend our results to chaos expansions in terms of the Poisson random 
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measu re, with the use of the relationship between the two chaos expansions derived bv lBenth et al. 
([2003). 



4 The Chaos Expansion with respect to the Poisson random measure 



ltd (|1956l ) proved a chaos expansion for general Levy processes in terms of multiple integrals 



with respect to the compensated Poisson ra ndom measur e. Note that it is trivial to covert the 
representation to iterated integrals as done bv lLokkal ([20041 ). The compensated Poisson measure is 
defined to be N (dt, dx) = N (dt, dx) — v (dx) dt, where v (dx) is the Levy measure of the underlying 
Levy process, X, and 

N(B) = #{t: (t,AX t ) B G B ([0, T] x R ) , 

where B ([0, T] x IRq) is the Borel a -algebra of [0, T] x IRq and Mo = {0}, is the jump measure 
of the process and its compensator is known to be 



E [N (dt, dx)} = v (dx) dt. 



Let / be a real function on ([0, T] 
(h,x{) , (t n ,x n ) is defined by 



Its symmetrization f with respect to the variables 



/ (tl, Xl, t n , X n ) I / , / (^7Tl 5 I •••) t-K n ) ^TTn) ) 

n! — ' 



(22) 



where the sum is taken over all permutations it of {1, ...,n} . f is said to be symmetric if / = /. 
The definition of symmetrization is used to represent the CRP in terms of multiple integrals instead 
of iterative integrals. 

4.1 Pure jump case 



We first consider the representation for pure jump Levy processes as in lLokkal (|2004l ). Let 
Li ((A x v) n ) be the space of all square integrable symmetric functions on ([0, T] x R) n . In an 
iterative integral such as Q, the time variables t\, ...,i n are monotonia For ease of notation, we 
let 

G n = {(h,xi, ...,t n ,x n ) : < t < ■■■ < t n < T;xi 6 R,i = 1, ...,n} , (23) 
and let Li (G n ) be the space of functions g such that 



9\\L 2 (G n ) = 9 (ti,xi,...,t n ,x n )dt 1 v(dx 1 )---dt n i'(dx n ) < oo. 



For / e L 2 (G n ) , let 



Jn(f) 



f (h,xi, ...,t n ,x n ) N (dti, dxi) • • • N (dt n , dx r 
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For / G L 2 ((A x v) n ) , let 



InU)= / )N® n (dt,dx) = n!J n (/). 

./([0,T]xR) n 

Chaos Expansion for pure jump Levy processes: Let F be a square in tegrable random variable 
adapted to the underlying pure jump Levy process, X t . Iltd ()1956l ) proved that there exists 
a unique sequence {/n}^Lo wnere fn £ L 2 ([0, T] x M) n such that 



F = E(F) + £/ n (/ n ). 



(24) 



ra=l 



Benth e£ a?. 



(|2003l ) derived relations between the two chaos expansions, that is, between the 



expansion in terms of compensat e d pow er jump processes and the expansion in terms of the Poisson 
random measure. iBenth et al. I (|2003l ) showed that the compensated power jump process defined 



in satisfies the equation 

(t) = 

Hence, the CRP can be written as 



x*AT(ds,dx), 0<t <T, i = l,2, 



f=e(f)+y: e 

j=l ii,...,ij>l 

00 

= ^) + E E 

J=l u,...,ij>l 



00 ftl — 

() JO 

00 /• /<tl 

JR •/ 



/ ( , li ... jij . ) (t 1) ...,t,Od^ ) ...dy t f ) dy// l) (26) 



jr 



7 2 £ 1 

j x 2 x l 



x /(i 1 ,.„,y)(ti, tj)JV (dtj, dxj) ■ • • JV (dt 2 , dx 2 ) iV (dti, dxi) 

• / / ffj (tl,^!,...,^,^) 

JO il 

x N (dtj,dxj) ---N (dt 2 , dx 2 ) N(dti,dx{) 



E(F) + E 
i=i 



00 <>oo /■ /•ti- 

Jr Jo 



= E(F) + E Jj (Si) = E{F) + E n\J, (g 3 ) = E{F) + E h (&) , 
j=i j=i j=i 

where c/j is the symmetrization (defined in (|22|)) of the function gj given by 



(25) 



(27) 



9j (^1)^1) ■•■,tj,Xj) 

" I 



■<V ■ ■ ■ f{i u . ..,ij)(h, -,tj), on G r 



on ([0,T] xK)"-G„. 



(28) 



Therefore, by uniqueness, {/ n }^L in (|2"3]l is given by 
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This equation provides a simple relationship between the two expansions. From Theorem [3l of 
course, 

(X t+t0 -X t0 T= £ U^ t> A, t ,to + ■ (29) 

We can now use this relationship to derive a form for g n in terms of I n , H^ ta and c[ n } . Let 
%-l,s = |(^ij---!^) \ij S {1)2, ...,s} and^ =1 ij = sj . Since the length of a tuple must not be 
greater than the sum of all the elements in the tuple (because an element must be at least 1), 

n s 

I < s. By definition, we have Z n = [^J [^J /Q jS . So we can write 



=1 t=i 



1=1 s=l 8 n £Ki s 



where 6 n is the tuple yii", 1 with I elements which sum up to s. Therefore, we deduce that 

for F = (X t+to - X ti) ) n in (J2BJ), /(^...^(ti, ...,tj) is given by 

f (il> ..., ij) (t 1 ,...,t j ) = U^ tta . (30) 

By (|28p . we have then proved the following proposition. 



Proposition 6 For any pure jump Levy process Xt satisfying condition ([2]), 

n 

(x t+t0 -x to r = Y,iiU n) )+cS, 



i=i 

(n) (n) 

where g) is the symmetrization of the function g) defined by 

gl n) (ti,xi, ...,t h xi) 

Y,s=i X^„e/q s x i "' x j n 0„,t,cr' on G r , 



on ([0,T]xR) n -G n , 



where C^" and LTg 1 ^ a are defined in Definition (TJ 
The following proposition gives a more straightforward representation. 



Proposition 7 For any pure jump Levy process Xt satisfying condition ([2]), 

roo r rti- r ftj-i~ f fin e n e n 

(x t+to -x to r= £ / / / / •• / H ■■■< 2 ^ 

a cT jo Jr Jo Jm. Jo Jr 
x ng t)CT iV (dtj, d Xj ) • • • N (dt 2 , dx 2 ) N (dii, dxi) + Cg, (31) 
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where c[ n } and a are denned in Definition Q3 

Proof. This follows directly by replacing j {i Xr , n i j ){t\t ■■■,tj) in (f27|) by ff30|) . 



□ 



Note that both chaos expansions, that is, the expansion in terms of compensated power jump 
processes and the expansion in terms of random measure, depend on I n , H^ ta and c[ n ^ . From 
(|25|) . we note the relationship between (t) and N (ds, dx) . Because of the simple form of this 
relationship, we can use Theorem Q] to derive the explicit representation of (|3ip . 

4.2 General case 



We shall now discuss the general relationship betwe en the two represen tations. Iltol (j 19561 ) 



proved the chaos expansion for general Levy functionals. iBenth et al. I (|2003l ) gave the relationship 
between the chaos expansions in the case with both a continuous (Wiener process) component and 
a pure jump (Poisson random measure) component. In this general case, the stochastic integrals 
are in terms of bo th Brownian motion , W, and the compensated Poisson measure, N (•, •). Hence, 
to unify notation, 



Benth et al. 



42003 ) defined: 



Ui 



u 2 



Ui = [0, T] and U 2 = [0, T]xK 
dQi(-) = dW r (-) and Q 2 (.)=N(-,-) 
d (Qi) = dA and d (Q 2 ) = dX x du 

9 Ql (du (1) ) = J g{s)W {ds) and 

<7(« {2) )Q2 (du^) = C I g{s,x)N{ds,dx) 



Benth et al. 



where A (dt) = dt is the Lebesgue measure on [0, T] . Let F be a square integrable random variable 
adapted to the filtration generated by the underlying Levy process, X. 
that 

oo 

F = E[F]+J2 E J ™ *° 

n=l j 1 ,...,j„=l,2 



(|2003h proved 
(32) 



for a unique sequence g^ 1 '"'^ (ji,...,j n = 1,2; n = 1,2,...) of deterministic functions in the 
corresponding L 2 -space, L 2 (G n ) = L 2 (G n , ®" =1 d (Q^)) , where 

G n = { (u[ jl \ e Tl^Uj, : < h < ■ ■ ■ < t n < r} 

with = t if ji = 1, and u^^ = (t, x) if ji = 2, and 

Jn ('/;/ " : ) 



l Gn Uf\...^)Q n {duf l} ) ■■■Q jn (d«g»>) . 



,0'i) 
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Similar to the pure jump case, we can derive the explicit formula for the chaos expansion with 
respect to the Poisson random measure of a general Levy process, i.e. a 7^ 0. In this case, we have 



t t-t 



Y^(t) = a / dW(ds)+ / / xN(ds,dx) 
Jo Jo 



Y« (t)= / x l N (ds, Ax) , < t < T, i = 2, 3, .... 
Jo Jr 

To derive the relation between the two chaos expansions, we introduce the following notation: 

flW (ds, dac) = adW (ds) + f xN(ds,dx) 

Jr 



i? W (ds,dx) = / x*N (ds, ds) , i = 2, 3, .... 
Jr 

Hence the CRP with respect to the power jump processes can be written as 

F = E(F)+jr, £ nr... f : , Ali / / M>v;-.,n;;/'«ny: ' 



w+E E / / •••/ 



/(ii,...,^)^!'-'^') 



x i?^) (dtj,dxj-) ...i? (i2) (di 2 ,dx)i? (il) (dti,ds) 
From Theorem [3J 

(X t+t() - X to ) n = £ n e™!t,<7 5 6» n ,t,to + ^l," 



oo rti— rtj—i 



r(») 



= £ / / -/ " .,, „ 

fl ~r JO JO JO 

x (dfyday) ...i?^") {dt 2 ,dx) R^ n) (dti,dx) + C t ( ™ } . 
We have then proved the following proposition. 

Proposition 8 For any Levy process Xt satisfying condition @, 



(x t+to -x to r= £ / / .../ " l>) 



Jo Jo 



rr 

9 n ,t,a 



R^ n) (dtj^xj) ...R^ n) (dt 2 , da) i? (i i n) (dti, dx) + C t {n) 



x 



^( n ) ttOO 



where C+J and FT), , are defined in Definition [H 
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5 The explicit chaos expansions for a common kind of Levy functionals 

Note that we have only found the explicit representations for powers of increments of Levy 
processes. In this section, we explain how the explicit formulae for a common kind of Levy 
functionals might be obtained using multivariate Taylor expansion. 

Assume that a real function g, possessing derivatives of all orders, is such that 



F — 9 { x h,X t2 - X tl ,...,X tn - Xt n _j) , 



(33) 



where the indices < t\ < t 2 < ■ ■ ■ < t n are known and n is finite. By expressing F in terms of 
power of increments of X, we can use our explicit formula to obtain the CRP of F. For example, in 
financial applications, g corresponds to all pricing functions of contingent claims which depend on 
the underlying asset at a finite number of time points. Suppose {Xt , < t < T} is the background 
driving Levy process and time is now t = t n . Suppose the un derlying asset, {St, < t < T}, is 
given by the exponential-Levy model, see lCont Tankovl (|2003l . Chapter 8.4): 



S t = Soexp (X t ) , 



where Sq is the initial value of the underlying asset at time t = 0. Then, for example, we can 
represent F as the pricing functions of a number of contingent claims listed in Table 1. 



In ([33]), let xi = X tl , x 2 = X t2 - X tl , x n = X tn - X tn _ x . 
powers of Xi, we need to use the multivariate Taylor series, see 
the points Xi = 0,i = 1, ...,n to obtain such a representation: 



f q is not a linear combin ation of 



Jeffreys k, Jeffreys! (| 19881 ). about 



oo 1 

g(x 1 ,...,x n ) = J2 < j 



j=0 



n d 

^ Xk dxl 



g (x^, x' n ) 



(34) 



xi=0,...,<=0 



Note that this representation exists when g possesses derivatives of all orders at zero. To show 
typical elements in this representation, we note the special case of n = 2: 



°° J 1 



d 



d 



Xl ^~r + X2 Trr 

1 " ox n 



+ 



j=0 

5 (0,0) + 
1 



dx[ 



9 (4,4) 



xi 



dg 
dx[ 



si=0,3;' 2 =0 



2! 



2 d 2 9 



1 dxf 



4=0,4=0 dx[dx' 2 




+ x 2 2 



d 2 g 



dx' 2 



4=0,4=0 
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Name 


Formula 


Forward and future contracts on a 
security providing no income 


F t = S t exp (r (T - t)) = S exp (X t + r (T — t)) , 
where r is the risk free interest rate and T is the 
maturity of the contract. 


Forward and future contracts on a 
security providing a known cash 
income 


F t = (S t -I)exp (r (T - t)) 

= (So exp (X t )-I)exp (r(T-t)), 
where I is the present value of the perfectly 
predictable income on S. 


Forward and future contracts on a 
foreign currency 


F t = S t exp((r-r f )(T-t)) 

= So exp (X t + (r — rj) (T — t)) , where ry is the 
risk free interest rate of the foreign currency. 


Forward and future contracts on 
commodity 


F t = (S t + U) exp(r(T-t)) 

= (S exp (X t ) + U) exp(r(T-t)), 
where U is the present value of all storage costs. 


European call options 


F (t, S t ) = exp (-r (T - t)) E Q [(S T - K) + \T t ] , 
where K is the strike, T is the maturity, Q is the 
risk neutral measure and Tt is the filtration of S. 


Up-and-out barrier call options 


F(t,S t )=exp(-r(T-t))E Q 

where H is the barrier and 
M t s = sup {S u ; 0<u<t}, < 


(S T - K) + l{ M s <H } 
t < T. 


1 


Up-and-in barrier call options 


F(t,S t )=ex V (-r(T-t))E Q 


(S T -K) + lj M |>#j 




Down-and-out barrier call options 


F(t,S t )=ex V (-r(T-t))E Q 
where mf = inf {S u ; < u < t} 


{St - K) + l{ m s >H j 
, < t < T. 


; 


Down-and-in barrier call options 


F(t,S t )=ex V (-r(T-t))E Q 


(S T -K) + l{ TO s<,H-j 




Lookback options with a floating strike 


F(t,S t )=ex V (-r(T-t))E Q 


M| - St] . 


Lookback options with a fixed strike 


F(t,S t )=ex V (-r(T-t))E Q 


(M§,-K) + . 


Asian call options 


F (t, S t ) = «P(-r(T-t)) Eq (ELi ^ _ nR)+ ^ 



Table 1: The contingent claims and their pricing formulae to which Taylor's expansion can be 
applied at some values of St- 



Let (°) 

assume that 



I 9 9 

/! &e'. dx'. -dx'. 

31 J2 Ji 



. As in 



Corcuera al. 



x' x =0,...,x' n =0 



E E ^* *(«) 

'= 2 iiv-JiG{l,...,n} 



i?' < OO, 



(|2005l . Lemma 2), we 



(35) 



for all i? > 0. The multivariate Taylor series (|34p expresses i 7 in terms of sum of products of 
powers of increments of X. From Theorem [H we can substitute Xi, i = 1,2,. ..with the iterated 
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integrals with respect to the orthogonal martingales. 
For all F G L 2 (f2, F) having the form (|55j). then 



j=0 



1 







.k=l 



n i n 

= g (0, 0, 0) + J2 " *%-x) ^ (0) + ^ E " ^-0 2 ^ (°) 
i=i ' j=i 

' 31=132=1 
1 n 

+ 3!E(^-^-0 3 41(°) 



ra n 



31=1 32=1 
n n n 



+ 37 X Yl Yl Hh^h^a} ( X t n ~ X t n -i) { X h 2 ~ X t n -i 



jl=l }2=l 33=1 



Xt n ~ X * J3 -i J 9hfrjs (0) + • • • , 



(36) 



where (Xj. — Xt i _ 1 ) n 's are given by Theorem U] and we assume = 0. The sums converge for 
every u> G f2 because of ([35]) , 



Since < < < • ■ ■ < t n , the product of two iterated integrals with non-overlapping limits 
results in an iterated integral: if i < j — 1, u, v G {1,2,3,...} and (f>i,(f>j are the predictable 
integrands, 

J0 JO 



J0 J0 



since r\ > tj—i > ij, giving an iterated integral. Hence, we get a chaos expansion of F in terms 
of iterated integrals with respect to orthogonalized compensated power jump processes. 

Note that in some applications, it is only necessary to apply Taylor's Theorem directly to F to 
obtain a PRP representation. While the approach given in this section gives the CRP of F, each 
(Xf . — Xt- -,) m c onsist s of an infinite sum and therefore (|36|) is composed of two levels of infinite 



sums. 



Yip et al. I (|2007l ) applied Taylor's Theorem directly to obtain the PRP of European and 



exotic option prices for hedging. 



STATISTICS SECTION TECHNICAL REPORT TR-07-04 



24 



6 Simulations using the explicit formula 



To verify the theoretical results given in Section [21 we simulate the underlying Levy processes and 
compare the values of (Xt+t ~ Xt ) n with the value given by its chaos expansion. In simulations 
we apply the stochastic Euler scheme for the stochastic differential equations (SDEs) of general 
Levy processes. The rate of convergence of this scheme for Levy processes was discussed by 



Protter Talavl (119971). For an up-t o -date int r oduct i on to numerica. 



solutions of SDEs, see 



for ex ample Iffigham & Kloedenl (12002 ). Iffighaml (j200ll ). iKloedenl (j200J) and 
(11999J). 



Kloeden Platen 



The processes considered are Gamma process and a combination of Wiener and Gamma pro- 
cesses. For illustration, we run simulations for k = 4 and k = 9 in the pure jump case and k = 5 
and k = 8 for the combined case. The plots produced are shown in Figures 1, 3, 5, 7 in Appendix 
lEl respectively. In the second and fourth simulations, we set to = 0.0099 and to = 0.0019 respec- 
tively. These simulations substantiate our explicit formula for the CRP for to > 0. We see that 
processes generated using the CRP and those generated directly from the Gamma process jump at 
the same time points. The differences between the two are plotted in Figures 2, 4, 6, 8 accordingly. 
Note that the axis of Figures 2, 4, 6, 8 are in much smaller scales than those in Figures 1, 3, 5, 7. 
We deduce that the difference is due to approximation errors of the stochastic Euler scheme. The 
errors decrease with the step size A. In each of the Figures 1, 3, 5, 7, independent realizations of 
the Gamma and Wiener processes are used. 

7 Conclusion 

Levy processes were introduced in mathematical finance to improve the performance of some of 
the financial models which are based on using Brownian motion as the underlying process and to 
model stylized features observed in financial processes. The derivation of an explicit formula fo r 



the CRP has 



Benth et al. 



jeen o f the fo c us of consi derable study, f o r prev ious work, see 



Leon et al. 



(2002), 



(|2003l ). lLokkal (|200^ ) and lEddahbi et al. I (|2005l ). In this paper, we have derived 
a computational explicit formula for the construction of CRP of the powers of increments of 
Levy processes in terms of orth ogonal comp e nsate d power jump processes and its CRP in terms 
of Poisson random measures. iJamshidianl (|2005l ) extended the CRP in terms of power jump 



processes to a large class of se mimartingales and we have shown that in the Levy case, our formula 
complements the one given by IJamshidianl (|2005l ). Our explicit formula shows that the integrands 
of the stochastic integrals in the CRP of the powers of increments of Levy processes do not depend 
on the integrating variables nor the starting time, which makes the construction and simulation of 
the CRP much easier. The coefficients of the CRP depend on m^'s which represent the moments of 
the process with respect to its Levy measure. In this paper, we consider only Levy processes and 
their compensators are always of the form rriit. Using the same calculation, it is trivial to extend 
the representation to semimartingales which stochastic compensators have known representations. 
The CRP of the pricing functions for some common financial derivatives can be found by expressing 
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the pricing functions in terms of powers of increments of the underlying Levy process using Taylor's 
expansion. 



APPENDICES 



A A note on the Nualart and Schoutens representation 



Nualart fe Schouteng ( 2000 ) derived the basic result for representing (X t +t„ — X to ) k when t = and 



k = 2. In the proof of the CRP. lNualart fc Schoutens! (|2000r ) made use of Proposition 2 in their paper, given 
i n Sec tion l2.1l and the following equation derived from the Ito formula (equation (5) in 

daood)): 



Nualart fe Schoutens 



(Xt+to ~ X tQ ) 



k{k-l)[ {X t+t0 -X t0 ) k 2 t- I a d(X. 



L s+t 



J2 (*) I ° (X^-X t0 ) k ^dY^+^(^rn 3 t(X t+t0 ~X t0 ) k ^ 

k - 1 /U\ ft+to 

V . )m J / s d{X s - X to ) k ~ 3 + m k t. 

• =1 V7/ J tg 



(A.l) 
(A.2) 

(A.3) 



There is a small inaccuracy in this equation and we provide the corrected one necessary for the derivation 
of the explicit formula. The second term in (|A.2[) should be 



.7 = 1 V ^ 



rather than X^=i (j) m jt i^t+to ~ -^io) 3 ■ The error propagates from equation (4) in 
pOOOf ). By integration by parts, J t * + *° (Xs- ~ X to ) k ~ j ds should give 



Nualart fc Schoutens 



k-l 



j'=i 

rather than the term 
fe-i 



£( • )m s (t + t )(X t+tQ -X to ) k - i - Y, 
J/ i=i 



t+to 



sd(X s ~X t0 ) k 3 +m k t 



k-l 



Y( )m jt (X t+t0 -X t0 ) k -i-Y 
j=i U/ i=i 



t+t 



sd(X s -X t0 ) k j +m k t 



Nualart fc Schoutens! (|2000l . p. 114). Omitting to makes the constant term of the representation 



stated in 

not equal to the expectation of (X t +t — X to ) k since it depends on to. Equation (5) in 



Nualart fc Schoutens 
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(|2000h should in fact be: 



k(k-l)[ (X t+t0 - X to f- 2 t - I sd (X s+t0 - X t0 ) 



E (*) / t+t " - + fj Qm, (t + t ) (A t+t0 - A 



3 

fc-1 

E 



t + t 



sd(X s -X t0 ) k 3 +m k t 



(A.4) 
(A.5) 

(A.6) 



As an illustration of this representation, we derive (G t +t — G to ) 2 using (| A. 1[) - (|A.3[) to inspect the constant 



terms. Note that Gt is a Levy process with a 2 = 0, so the terms in (|A.1[) are equal to zero. We have 

rt+to rt+to 

(G t+t0 ~G ta f = 2 (G tl - - G to ) dG^ + / d&fj + 2 mi t (G t+to - G to ) 

J t J t 



rt+to 

2mi / h d(G fl -G to ) + m 2 t 

Jtn 



to 
t+to 



Gill - G^l + mi {h - h) dG\:> + I d<?£ 

J Jto 
t+t 

(1 



, t+t 

M 1 ) , / AM 2 ) 



2mi / t\ d 

'to 



rt+to rt\ — rt+ta 

= 21 I dG^dG^ + / dGf + 2mi (t - t 



to Jto 

+ m 2 t 2 + m^t — 2m 2 tto- 



to 



G\\> + rmh 

t+t 
to 



m 2 t 



The expectation of 2 J t c J t 1 dGj 2 dG^ + J t "dG^ + 2mi (t — to) f tg "dG^ is zero since the compen- 
sated processes G| 1} and G {2) have zero means. We see that m\t 2 + m 2 t — 2m 2 tto depends on to which in 
fact cannot be the expectation of (Gt+t ~ Gt ) 2 since the increments of G* are stationary. Starting from 
([Of - (rA~6|) , we can find that 



rt+to rti — rt+to rt+to 

(Gt+to ~ Gto) 2 - 2 / / dG«dG« + 2mrt / dG« + / ^ ■ ^ 2 

J to J to J to J to 



dG^' + m\t 2 + mat. (A.7) 



B Proof of Proposition [2] 

We prove this result using strong induction. Clearly, the proposition is true for k = 1 and 2. Assume 
the proposition is true for k — n, where n is an integer > 1. Then for k = n + 1, firstly we prove that the 
sum of the indices of all the m 9 's appear in each of the terms of Cj (given by Proposition [T]) are equal 
to n + 1. We have 



Gr^E^D-^r^-E 



n + 1 



h dG t ( " +1 J *> + m n+1 t. 



(B.l) 



By the induction step, the tuples of the indices of all the m g 's appearing in each of the terms of G\ 



(n+l-i) 



are elements of C n +\-j defined in (|12p . Since we have irijC\ 



(n+l-i) 



appearing in the first term of (jB.l 



mj-G tl 



(n+l-j) 



in the second term and m n +i in the last term, it is clear that the tuples of the indices of all 
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the m g 's appearing in each of the terms of c\ n+1 ' are elements of C n+ \. Now from (|B.1[) , the first term 
can be proven to be 

'n+l\ ,^( n+ i_j) 



E 

j=i 



.1 



ijtC 



E 



^= 1 0„ +1 - i =(4" +1 -",4" +1 -".-4" +I ^ ) )e^ . 



E 



l\ V 1 



) *2 



(j'r.+X-j J> n +l-j ] , 

I Pi )P2 ' — )P„+l-j 



n 



ge0„ + i_jU{j'} 



/Z+l 



and the second term can be shown to be 



E 

i=i 



n + 1 

i 



m, / ii dC 



(n+l-j) 



E 



E 



Z! V* 1 



,«2 



-.7 



x I Pi >P2 ) -iPn+l- ' 



n 



iz+1 



Z + l 



Hence, 



(n+l) 



E 



E 



^= i ^ +w =rii" +i -«,4 n+i - ,) ,...,i, (n+i - j) )6£n +w 



1 An+l-j) An+l-j) A, 

> *2 ' — ) *j ' >/ 



X „^+i-i } ! 

x \Pl 1P2 > —>P n +l-j 



n 



/??,, 



g£0n+l-jU{j} 



Z + l 



m n+ it. 



Next we are going to prove that 

E 



(1+1) 



1 17™+!) ,-("+!) ,"("+!) \| 



n 



iZ+1 



E 



E 



(Z + l) 



\ (.(n+l-j) .(n+l-j) An+l-j) A , 



(Pi .Pa .->P„+i-j 



n 



9£*n+l-jU{j} 



+ m n+ i*. 



(B.2) 



On the R.H.S., we are adding a j to each tuple 



such that 4" +1 J) + 



j = n+l. Suppose <^ n +i = [i\ ,i 



= (M+l) •(«+!) An+1)\ 



pie fijf 



'2 1 ' 



J has one extra element compared to the tu- 



(n+l-j) An+l-j) An+l-j) 
, t 2 j •••) 



Z+l -(n+l) 



and otherwise they are the same. Since ^2 q=1 iq 



= n + 1, to 



STATISTICS SECTION TECHNICAL REPORT TR-07-04 



28 



obtain ( i^ 1 ^ 1 ', i^ 1 ^ 1 ', ) fr° m ? ^4 n+1 ^ > ■••> "''')' we are adding an element j to 



.(n+l) ,.(ti+1) -(»+i)A 

... y 

the latter such that the sum of the tuple is equal to n+l. Suppose there are r distinct value (s) 

+1 \^ +1 \-Jr + X 1] ) and let 
Note that Y^ q =i fq 1S equal to 



fAn+l) .(n+l) -(n+l) 
1*1 '2 !'"! 4 i+l 



Let ari, X2, x r be the distinct values in (i^f 1 



+1) ,-("+!) ,-("+!) 
, i 2 , t J+1 



/i, i = 1, r be the number of times Xi appears in (i^ 
the length of the tuple , , ...,4+T^) . that is ' Ej=i /, = I + 1- Since 4" +1) , ...z^) 

can be obtained by adding an element j to a tuple i 2 n+1_3 \ i[ n+1 3 M whose elements add up 

to n + 1 — j, j can take one of the r distinct value(s): x\,X2, ...,x r . For example, suppose j = Xi, then the 
corresponding term on the R.H.S. of (|B.2[) is 



(n + l)! 



(an!/ 1 (x 2 l) h ■ ■ ■ (xA f- 1 ■ ■ ■ (x r l) fr x^l •••■:/, 1 ;! • • • f r \ 
Summing up all the possible j s {x%, X2, ay} , 



n 



J+l 



Eg=l / q 



E n/i 



(n+l)! 



^ (xx!) J1 (a; 2 !) /2 • ■ ■ (x^ 1 • • ■ (ay!)'' ^!/i! • • • (/< - 1)! • • • / r ! 

(n + l)! 1 
(a;i!) /l ---(a: r !) / '/i!---/r!EU/9 

(n + l)! 



n 



d+i 



n 

<?60„ + l 



* i+1 E/* 



i=l 



(a:i!r---(a;r!r/i!---/r! 



n 



f ;+i 



For the case n +i = ^ n+1 ^ , it is clear that the L.H.S. of (|B.2[) is equal to m n+ it. Hence, by applying the 
same argument to each possible tuple 4+1^) e Ai+i, we have proven (|B.2[) and therefore 



C, 



(n+l) 



E 



+1) ,-("+!) -(n+l) \ j 
, t 2 j ■ ■ ■ ; 



n 



C Proof of Theorem [T] 

We prove the result using strong induction. Firstly we need to consider defined by equation (TIT 
We need to know what tuples are in Tk+i but not in X^, and these correspond to those elements adding up 
exactly to k + 1. Let Jk+i be the collection of these tuples, that is, Ik+i =Zk U Jfc+i- We have 



Jfc+i = |(«i,«2, I j £ {1,2, fc + 1} , i p e {1, 2, ...,fe+ 1} and = fc + 1 j . 

To construct J7fc+i from we can simply add an element to the end of each tuple in Ik so that the elements 
of each new tuple add up exactly to k + 1, and finally including the tuple (fc + 1) in J^k+i- 

We are going to prove by strong induction that (G t+to - G to ) fc = Ee fc ei fc ^6%^0 k ,t,t + for any 
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non-negative integer k. For k = 0, clearly both sides equal 1. For k = 1 and 2, it can be checked 
easily that the proposition is true. Assume the proposition is true for k — Q, 1,2, ...,n, where n is a 
positive integer. Note that it is sufficient to prove the representation for G™ +1 only since we can always 
let Gt+t — Gt a = F t , which is also a Levy process and we have, the z-th power jump process of {Ft, t > 0} , 
F t {t) = - G^j for i = 1,2,3, .... Since both {F t ,t > 0} and {G t ,t > 0} are created by the same 

infinitely divisible distribution, the compensators for their i-th power jump processes are both equal to rait. 
Hence, we have the i-th compensated power jump process of {F t , t > 0} , 

F«=G« o -G^. (C.l) 

For k = n + 1, by (|A.4|) - (|A.6|I . 

or 1 - E ) / <%t--' + E { ] ) miter 1 -* 

Firstly, we want to prove that all the stochastic integrals in G™ +1 is of the form 5e n+ll t,0) where n +i S 1n+\- 
From (|C.2p . it is clear that the first term is the only term introducing new stochastic integrals which are 
not in T n . The general term of the stochastic integrals in the first term is 

'g^-UG^, i = l,2,...,n + l. (C.3) 

o 

By assumption, 

G ti- - 2^ 11 n+ i-,',ti 6e «+i-i.*i.° + C *i , J - l,2,...,n+ 1. 

9„+i-j£l n +i-j 

When j = 1 in ([03]) . we have /„' G^_dG^ ; 3 meaning that we are adding alto the end of all tuples in I n . 
Since by definition 

J„ = |(ii,i 2 ,...,ij) | j S {l,2,...,n}, i p e {l,2,...,n} and E i p < n| , 

we know that the sum of the elements of the new tuples we get from adding a 1 to the end of each tuple 
of T n is less than or equal to n + 1. Similarly, when j = 2, we have J* G™ i ^ 1 dG^' ) , meaning that we are 
adding a 2 to the end of all tuples in X„_i and since by definition 




we know that the sum of the elements of the new tuples we get from adding a 2 to the end of each tuple of 
X n _i is less than or equal to n + 1. We can continue the same argument until j — n. When j = n + 1, we 
have J *dGj™ +1 ''. Since T n D T n —\ 3 •■• D I2 ^> 2Ti , the above way of introducing new stochastic integrals is 
the same as adding an element to the end of each tuple in I n so that the elements of each new tuple add 
up exactly to n + 1. Hence all the elements in J n+ i have been created and since U Jn+i, w e 
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have proved that all the stochastic integrals in G™ have the form Sg n+1 t,o, where 9 n+ i € 1 n +i 



By definition, G t is the term in G" +1 not containing any stochastic integral. Hence it is correct to 
write c[ n+1 ^ as the final term. 

Finally, we want to consider the coefficients of the stochastic integrals, that is, we are going to prove 
Proposition [3] By assumption of the induction step, we have 



or 1 = E 

n 

+E 

3=1 

71 

-E 

3=1 
n 

+E 

3=1 
n 

-E 

3=1 



n + 1 

j 

n + 1 

j 

n + 1 
j 

n + 1 
j 

n + 1 
j 



E 

, ex 



+i-jti„+i-j 



771 jt 



E 

e„ + i_j-ex„ + i_ 

E 

971+1-36X71+1-; 



1 V+i- J ,ti 6e »+i-^*i.° dG *i 



ii e„+i-,-,t c ' e «+i-i.*. 



ti d 



li e„+i- J ,ii^ e ™+i-3' t i.° 



'1 ^E 

3=1 



n+1 
J 



mjtC} 



(n+l-j) 



m, ti d 



C, 



(n+l-i) 
ti 



dG^ +1) +m„ +1 t 



L x + £ 2 + is + i 4 + i5 + ie+ / dG|™ +1) + m n+1 t 



(C.4) 



Let /Q. s = • h) \ij £ {1) 2, s} and Y] l j—i ij — s\ . Since the length of a tuple must not be greater 
than the sum of all the elements in the tuple (because an element must be at least 1), I < s. By definition, 

n s 

we have 2 n = [J [J 7Cj jS . For any 9^ s e fCi tS , let 6i jS = 



i 2 »•••)'; 



It is obvious from Proposition 

8=1 1=1 

[I] that C[ k ' has the form c\ k ^ = q^ + q[ k h + t 2 + ... + q k t k . Note that q^ is non-zero only when 
fc = 0. When k — 0, by definition c\ = 1, so we have q^ = 1. We need to find out the recursive 
relationships between the gr 's. From (fTTj) . for fc > 1, 



(fe) + , (fe) + 2 , 
<7l * + Q2 1 + 



^* fe = E 

3=1 
fc-1 



<7i 



1 * + <?2 * H l"9fe-j * 



3 = 1 



e : k / *i d + + . • • + qtptx-i 



k-l 



3=1 



k-l 

E 

3=1 



1 ,NV + ?,N)^ 



a^--"f -I I ^ ff ( fc -3')+fe-i+i 



By comparing the coefficients of t, q[ k ^ = m k - By comparing the coefficients of t r , r = 2, fc, 
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To ease notation, we let 



i 

ti- rti-i- 



Jo 



h = I - I dGl- ■ -dG'- 1 , I 3 = / I 2 dG[, q, = 



Note that it is only for simplicity in writing out the equations. When doing calculation, we should always 
use the long but clear notation. So, we have 

n n+l—j s ,t n+l-j—s 

L * = E E E E v> E q-*nid^. 

j=i s =i ;=i e,, s e/c,, s Jo t«=o 

n n+l—j [ s s n+l—j—s 

L 2 = Y, m i E Wi-^E E f^ + ^-^E E F 2 E q-^ +1I 3 

i=l s=l [ fcl 1=1 8,,eK IjS w=l 

n ra+l-j s { [ l 1 If 1 

^3 = - E TO J E E E F 2 I qo / til 2 dG^ + - qi t 2 I 3 + - qi / i?I 2 dGl 

J= i s =i 1=1 1 Jo 2 2 Jo 



^=E( ) / E e^^d^ 

,•=1 V 3 J J ° w=l 
n , , 1 \ n+l-j 

i5=E t h E 



J — 1 V 7 U> — 1 



*» = -E( ! V* E 



j— 1 x 7 w— 1 



10 + 1 



Next, consider L\ and L 3 . Let u,v £ {1, 2, n — 1} and u + v < n. In ii, when j = u, s = v (hence 

s <n-j), 

v r t / _i_ 1 \ I 1 ( n+ l—u -v ~) 

^ = E E / ,fn + i u 1 m n +i- M -«ti + E } *r 

tl_ ^^.7 ^^ "d4;^ ) •••d4:" ) d^. (c. 6) 

Since s = v, I £ {1, 2, , we have by definition ( , i 2 ''" , ) G In i 3 , when j = n+l—u — v 

(hence j e {1, n — 1}), s = u + v (hence s = n + 1 — j) and i^ i,s = u (hence i^' s < s), 

\ - \ - (n+1)! 



«=i e, „+„€*:, u+ „ • • • ^ {n + l-u-v)\ 

x/' il /'-r.../"-d4' : ")..-a 6 p')a^'. 

Jo Jo Jo Jo 

Since s = u + v and i^ l,s = u, X)p=i*p ' l + " = w ' we nave by definition i^'" -1- ", ■■■,if l _li +v ^j S J^. 
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in Li and 



v- /■* (n + l)l 1 

Jo Jo Jo 



u+v , . 

E\ - (n + 1) 

-Oi.u+vrOt 



1=1 e, u+v zk, u+v h- ' " • • • (n + 1 - « - «)! 

Jo Jo Jo Jo 
cancel each other. So we now have 

n+l— j [ s n+1— j — s 



j=i [ s=i [ ;=i e l , s eK liS Jo w=2 



Since g fe) = for fc > 0, 

n+l-j [ s ( 1 1 /-t 



= - E "M ^^-n E i i^n-i} E E F M o^* 2 ^ + ,qi / *?i2 dG 

i=l I -=1 I I=10,,.€/C ( ,. 12 2 7 ° 

n+l-j-s r 1 ft 1 1 

+ E qJ^T^Ia + ^r / C +1 I 2 dGa 

^2 I «> + 1 W + 1 7 J J 



Next, consider L3 and £4. Let u <E {1, 2, n}. In L3, when j = n + 1 — u (hence j G {1, 

5), we h 

(n+l)! 

u! ( 

In L4, when j — u, we have 



(hence s = n + 1 — j) and z" ! s = u (hence i® l,s — s), we have 



L 3 = m n+ i- u \ 1 / *i d ^ } - 

u! (n + 1 - u)! j" 



L 4 = 



tu=2 



Hence the terms 

(n+l)! f\ ,Mu) 
m„+i_„— — — / h dG \ > 

in L 3 and 

rt 



V u ) J TO ™ +1 1 *i 
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cancel each other. In L 3l since the terms where (j — n) and (j < n — 1, s = n +■ 1 — j) get cancelled, we 
can sum j from 1 to n — 1 and sum s from 1 to n — j. 



n— 1 n—j s n+l—j—a 

^ = -E"^EE E F 2 E q 

3=1 s=l 1=1 6i,,eKi,. u>=l 
n-1 f / . , x ,( n+l-J -, rri+2-j-u; 



10 + 1 



1 



E 



+ 1 - j\ (n+l-]-z) 

z ) m zq y w -i 



t? dG« 



by ([CD. 

Consider L4 and L3. Let u G {1, 2, n — 1}, u G {1, 2, n — 1} , u + v < n, x G {1, 2, v} and 
hence x + u < n. In L4, when j = u, w = n + 2 — u — v (hence w G {2, ...,n + I — j}), z = x (hence 
ze {1, ...,n + 2- j -«;}), 



i4= 1 „ Ji n + 2-rx-A x ) m ^ + i-u- v h dG tl . 



(C.7) 



In L3, when j = x (hence j G {1, ...,n — 1}), s — u (hence s < n — j), i ; 1,s = u (hence j ; , s = s), w 
n + 1 — u — v (hence w < n + 1 — s — j), 

T _ (n + I )! 1 (n+l-u-x) 

^3 — m x -j—, — ; - - rrrSn+i-ti-u 



u\x\ (n + I — x — u)\ 



n + 1 — u — v 



n+2-u-v j ,^q( u ) 



1 



n + 2 — u — v J 11 n + 2 — u — vJq 1 

where the 2 nd term cancels (|C.7|) . So now we have 



n+2— 11— v ^q( u ) 



L 4 = 0. 



ri— 1 n—j 



(1) 

tl 



J=l 8=1 k 

s ^ n+1— j— s 

+ l{2<s<n-j} E E 1 1 (if'.«<8) F2 E * 

i=i e hs eKi, s I 1 J ««=i 

(n+1)! 1 
+ {#■'=.} (n + l-i-s)! 



U> + 1 W + 1 7o 



i™ +1 I 2 dG< 



Next, consider Li and L 3 . By the equation for q£ 3 s given in (|C.5p . we have 

n ( n+l-j ( s „( n+l-i-s _^ 

L i = E \ hj<n-i} E ^ 1 {«<n-i-i} E E / F2 E - 

i=l [ a=l [ !=1 O^eC,,/ ™=2 



n+2 — j— s— w 



x E 

z=l 

+ l{s=n+l-j"} 

Let u G {1, 2, n — 2} , de {1,2, ...,n- 2}, u+w < n-1, x G {1, 2, w} , /3 G {1,2, + 1 -x}. In L x , 



£ E /'FxIid^l + l^Cn + l)/' T ■ 
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when j — u, s — n — u — v (hence s G {1, n — 1 — j}), w = x+l (hence w G {2, n + 1 — j — s}), z = [3 
(hence z £ {1, ...,n + 2 — j — s — w}), 



^= E E 

f=i ej.n-a-.ex;. „_„_„ 



(n + 1)! 1 1 /u + 1 



o if-—" !#"™ !••• if 1 -— "!u!(« + + 1 V /3 
x^-V / ' ••• dci^ ■■■^i ,m) d6^. (C.8) 



o Jo 



By definition, since s = n — tt — w and £ G {1, 2, n — u — v} , fij l,n " "ji^'" " " "J G ^J n — u —v 

In £3, when j = (3 (hence j £ {1, n. — 2}), s = n — v (hence s G {2, n — j}), i® l,s — u (hence 
i® 1 ' 3 < s), w — x (hence w G {1, n + 1 — j — s}), 

™ ^ (n + 1)! 1 



l=i e,„_„e/Ci 



Jo Jo 



Jo 



The final term in L 3 

(n + 1)! 



m ^ E E ^..— ^..-.j . . . ^ r (« + !-/?)! 



x g (, + i-« _1_ / r i / • .../■• d( M> ^..-dG^- 1 ^ 



t /-tl- ft!-!- r,%n-v\ f.»l,n 







i + Uo 1 

clearly cancels (|C.8|) in Li. So now we can write 

1+1— j s „t 



ti • 



£1 = E 1 {j<»} E hs=n+i-j} EE/ FiIx d< ^ 

i=i s=i 1=1 9i, s eKM 



3=1 s=l 



j— s 



+ i{2<.<n- 3 -} E E iV' s <4 F2 5- ^^r+r <lu+l13 

z=i e t a eKi s (. 1 ' w=i 



(n + 1)' 1 



{*i'"'=»} s!j! (n + l-i-s)! 
We can now simplify it as 



n—l n—j s n+1— j — s 

E^EE E F 2 E 

3=1 s=l J=l 6i s elCi s iu=l 
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All together, we have 

* n+l-j 



£l = E 1 {3<n} \ E k^n+l-j} E E / Fill dG 



CO 

tl 



i=i ( «=i /=i Si,,eKi, 3 

n n+l—j I s s n+1— j— s 

^ = E m ; E ^{-n+i-iiE E ^ii 3 + i {s <„_, } ^ E p 2 E q^ +1I 3 

i=l s=l [ 1=1 9i >e eKi, B 1=1 8 l:S £K ls w=l 

n— 1 n—j s n+1— i— s 

^ = -E-.EE E f 2 e ^^ ,+1i 3- 

3=1 s=i 1=1 e ls eKi iS w=i 

L 4 = 0. 

n / , i\ n+1-3 



Since at the beginning of the proof, we have already showed that the stochastic integrals of G™ are of the 
form Se n+lt t,o where 6 n +i,t G 2~n+i- We are now going to show that the coefficient of each Sg n+1) t,o is ni n+1 ' . 

Consider • • • /*'" dG^ 1 ■■■ 1 dG^ where 9 1>S e /C,, s , j £ {l,2,..,n}, s = n + 1 - j. 

This stochastic integral only appears in L\. And its coefficient is — — l" +1 - >! e . And from ([Ml) , since 

n + 1 — s — j = n + 1 — (n + 1 — j) — j = 0, 



n (n+i) _ (n + 1)! ^(o) _ (n+jQ! 
u /.9(.a .«i. s a _ .8,...o,.. .9,, u * - 



,» 2 >— > l i >3) % x h 2 ! • • • i l !j!(J! i 1 h 2 ' ! • • • i ; !j! 

since Gj -* = 1 by definition (fTTj) . Hence we have proved that the coefficient is given by TL^t^ 

(*i ,s .v ■.' j 

Next, we change the summation sign of j and s in L%. 

A f^A x - (n + 1)! 



'■ s !^- s !---z; ! ' s !(n+l-s)! 



n—s s n+1— j— s 

+EE E ^ E 

j=i (=1 !iS e/c,, s tu=i 
Similarly, by changing the summation sign of j and w, we have 



s =i {i=i8i, e eK,, 3 t 1 l - s li 2 l - B l---i l i - s l{n+l- s)\ 

n — s n-\-l—w — s s 

"E E E E 

w=i j=i 1=1 «, ,eK| . 
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By m, ^ r;^~ w - s J ^:;l ) ^ +1 - s - j) = gft 1 -', so we h ave 



n | s 

L2 = \ E E m n+lst 

a=l [ 1=1 6 t , s eKi,. ' 1 



(n + l)\ 



n — s s 



EE E 

w=l 1=1 e t „e/Ci „ l l :t 2 



(n+1)! 



(t» + 1) 



(n + l-s)!^ 1 



<Wl £ *-3 



Changing J2Z=i to £™=2 S > we have 



Similarly, 



L2 = ^ \ E E m n+i- s t 
n+l— s s 



(n+l)! 



E E E 



1 -'2 

(n+l)! 



. ! (n + 1 - a)! 

w=2 1=1 e,, s G/C(, s n "2 ■ H ■ x ' 



•^ ! ' s !(n + l-.s)! 



n— 1 n — s n+l— to— s s 

L 3 = -EE E ^E E F 2 q.^ I ^i 3 . 

s =i tu=i j=i ( = i ej^e/Ci,, 



By^^jE^i 



+ 1— to— s (n+l— s)! (n+l — s— j) (n+l — s) 



3!(n+l-j-s)! m J^ 



, so we have 



n— 1 n+l— s s 

E E E E 



(n+l)! 



w — 1 



(n+l-s) ri 



s =i w =2 j=i e h3 eKi, a l i -H ■ H ■ K ' 
For s = 1, the stochastic integral J" Q * dG^ appears in both L 2 and i 3 . Its coefficient is given by 



J2 (n + 1) wq£h w + m n ( n + l)t-(n+l)J2(w-l) q^t v 



(n + l) 



= (n+l)Cj 



(n) 



(n+l) _ (n+l)! ^(n+1-1) = 



By (HI, 

L « ~ (n + 1-1)!"* 
For s £ {2, 3, n — 1}, the coefficients of the stochastic integral 



(n + 1) G, 



(n) 



JO 



dG 



ti 



•dG t v 2 MG^ 
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is given by 



m „ t (n+1)! + "v * + (^-i-.) t » 

z^k^ 3 !- ••«,'■■! (n + l-s)! w=2 • • • [n + l s)l 

_ (n+l)! (w-1) (n+1 _ s)w 

w=2 l l -H ■ l l ■ v ; 

E ( n + 1 )' 1 Jn+l-s) f w 
AV"I...A'I (n + l -a)! 9 " 

_ ( n + !) ! 1 -(rt+l-g) _ n (n+l) 



i?'-!i*"!...if l -!(n+l-«)! * (<!' 
by (|14p . For s — n, the stochastic integral appears in L2 only and its coefficient is given by 
mif + _ (» + !)» C (D _ n (n+D 



The stochastic integral J t dGj™ +1 ' > appears only once in G™ +1 and its coefficient is equal to one. By (fT4|) , 

(n+ i) _ (n+l)! (o) _ 
U (n+D- (n + l)! ' ~ 

Finally, we have to show that L 5 +L 6 +m n+1 t = c[ n+1) . By gUJ, ^ ("^K?™ +W) = 

n / , \ n+l— j n n+l— -10 / , 1 \ n 

L «=£( B t H S e +w) ^ +1 =E E n V,e +i - ) ^ +i -E(-+ i )^ + + i 1) ^ +l - 

j=i V J / w=1 w =i j=i \ J J w= x 

n n-\-l—w / 1 1 \ n 

^ = -E E nV^-^^^-E^ 1 ^ 1 - 

10=1 J = l 10=1 

Hence 

n n+l n+l 

L 5 + L 6 + m n+1 i = £ g^V* 1 + "Wi* = E + = E = C nl- 

w—l w=2 w — 1 

Thus, we have proved that 

As explained in (|C.1[) , since F f = Gt+t — Gt is also a Levy process, we can write 

^ +1 = EE E n K } [ C- r -1 "^^ ' s) 

s =i ;=i is £lC, , 



In . . . . ,1 , 



JO 
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and since cLF 1 /^ = d (g^^ — GtoJ — dG(+ to , by changing of variables, we have 



n+l s 



(Gt+to - G ta ) n+l - X! X! X! n e?T,t 



) 



^2 ^d^+i' 4 ^ + c < 



(n+l) 



Therefore, by the principle of strong induction, 



(Gt+t - G to ) k — ^2 ^-e k ,t^s k ,t,to + G, 



for all non-negative integers A;. 

D Proof of Proposition [4] 

We prove by induction. Assume the proposition is true for all k > n. Now, consider n 

n n ( I — 1 



yin+X) = H (n+X) _ £ an+i i y(l) = ff (n+X) _ £ + £ ^tfW 

n I n 



i=i fe=i 



fe=i 



which completes the proof. 
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Figure 1: G\ generated using CRP and directly Figure 2: The difference of the two series in Figui 



from the Gamma process. 



1. 
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Figure 5: A t 5 generated using CRP and directly Figure 6: The difference of the two series in Figure 
from the Wiener and Gamma processes. 5. 
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Figures 1-8: Solid line is generated using the CRP and the dotted line is generated by the Wiener and 
Gamma processes. Time step = 10 q 00 , a = 10, b = 20. In Figure 3, to = 0.0099; in Figure 5, o = 0.01; in 
Figure 7, t = 0.0019 and a = 0.02. 
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